function moment = hwozd_m( I,p,q )
    xs = (I.x) .^p;
    ys = (I.y) .^q;
    
    moment = sum(sum( (xs .* ys) .* I.image));
